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Properties of BaTiOs/BaZrOs ferroelectric superlattices with competing instabilities 
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Properties of (BaTi03)i/(BaZr03)„ ferroelectric superlattices (SLs) with n = 1-7 grown in the 
[001] direction are calculated from first principles within the density functional theory. It is revealed 
that the quasi-two-dimensional ferroelectricity occurs in these SLs in the barium titanate layers 
with a thickness of one unit cell; the polarization is oriented in the layer plane and weakly interacts 
with the polarization in neighboring layers. The ferroelectric ordering energy and the height of the 
barrier separating difi^erent orientational states of polarization in these SLs are sufficiently large 
to provide the formation of an array of independent polarized planes at 300 K. The effect of the 
structural instability on the properties of SLs is considered. It is shown that the ground state is 
a result of simultaneous condensation of the Lis polar phonon and phonons at the M point (for 
SLs with even period) or at the A point (for SLs with odd period); it is a polar structure with 
out-of-phase rotations of the octahedra in neighboring layers, in which highly polarized layers are 
spatially separated from the layers with strong rotations. The competition between the ferroelectric 
and structural instabilities in biaxially compressed SLs manifests itself in that the switching on of 
the octahedra rotations leads to an abrupt change of the polarization direction and can cause an 
improper ferroelectric phase transition to occur. It was shown that the experimentally observed z- 
component of polarization in the SLs can appear only as a result of the mechanical stress relaxation. 

PACS numbers: 63.22.-m, 68.65.Cd, 77.84.Cg, 81.05.Zx 
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I. INTRODUCTION 

Ferroelectric superlattices (SLs) represent a new class 
of artificial materials whose properties, such as the Curie 
temperature and spontaneous polarization, often exceed 
those of bulk ferroelectrics.— '^ An analysis of properties of 
ten ferroelectric SLs with the perovskite structure grown 
in the [001] direction have shown that they are ferro- 
electrics at low temperatures, no matter whether they 
are composed of two ferroelectrics, a ferroelectric and a 
paraelectric, or even two paraelectrics<^ One of the pos- 
sible applications of SLs can be non- volatile ferroelectric 
random access memory (FRAM) devices in electronics. 

The recent studies of the ferroelectric instability in 
free-standing (KNb03)i(KTa03)„ superlattices grown 
in the [001] direction showed that, with increasing n, 
the tendency to the ferroelectric ordering is retained in 
KNbOa layers with a thickness of one unit cell while the 
interaction energy between neighboring KNbOa layers 
exponentially decreases 4i^ It was shown that at n > 2 
the ground state of the SLs is an array of ferroelectri- 
cally polarized planes whose polarization is concentrated 
in the KNbOa layer, is parallel to the [110] direction, and 
weakly interacts with polarization in neighboring layers. 
The use of such arrays of quasi-two-dimensional ferro- 
electric planes as a medium for the three-dimensional 
data recording can make it possible to achieve a volume 
recording density of ~10^^ bit/cm^. 

Unfortunately, the implementation of the above- 
described ground state in the (KNb03)i(KTa03)„ SLs 
requires very low temperatures (~4 K) because of the low 
ferroelectric ordering energy and low height of the poten- 
tial barrier separating different orientational states of po- 
larization. That is why one of the aims of this work was 



to study the feasibility of similar quasi-two-dimensional 
structures and the ways to improve their characteristics 
in BaTiOs/BaZrOa ferroelectric superlattices with an ac- 
tive layer of barium titanate. 

In BaTi03/BaZr03 superlattices, the tension of the 
BaTi03 layers as a result of the epitaxial matching 
with the BaZrOa layers having a larger lattice parame- 
ter should stabilize the Amm2 ground state necessary to 
implement the quasi-two-dimensional structure* in the 
BaTi03 layers. However, since the simultaneous com- 
pression of the BaZr03 layers can cause the polarization 
normal to the interface to occur in them,— the search 
for the ground-state structure in these SLs requires a 
detailed analysis. Furthermore, the instability of the 
BaZrOa structure with respect to the octahedra rota- 
tion2r— can result in the competition of the ferroelectric 
and structural instabilities in the SLs, and this problem 
also calls for a detailed study. 

The properties of BaTi03/BaZr03 SLs have been 
studied previously both experimentallyiSrJ^ and theo- 
reticallyf2, The initial interest to these SLs was caused 
by giant dielectric constants observed in themJ^ Room- 
temperature observations of ferroelectric hysteresis loops 
in the SLsiiii^ii^ suggested that the ferroelectric tran- 
sition temperature is them exceeds 300 K. The spon- 
taneous polarization nonmonotonically varied with the 
period of SLs and depended on the electrode configura- 
tion (in short-period SLs with electrodes on both sides of 
the film, hysteresis loops were absent ^■^). These features 
have already been discussed in Ref. 3; and were explained 
by the rotation of the polarization to the layer plane. 
The recent study of Raman spectra and dielectric prop- 
erties of BaTi03/BaZr03 SLs revealed the appearance 
of the z-component of polarization in barium titanate 



layers in the SLs grown on MgO substrates. ^^ Interpre- 
tation of these results requires the study of the influence 
of substrate- induced strain on the groimd-state structure 
of the SLs. 



II. CALCULATION TECHNIQUE 

The (BaTi03)i(BaZr03)„ superlattices considered in 
this work are periodic structures grown in the [001] 
direction and consisting of the BaTiOa layer with a 
thickness of one unit cell and the BaZrOa layer with 
a thickness of n unit cells {I < n < 7). These SLs 
were modeled on supercells of IxlxL unit cells, where 
L = n + 1 is the SL period; in modeling the structures 
generated by instabilities at the M and A points at the 
boundary of the Brillouin zone, the primitive cell vol- 
ume was doubled. In this work, particular attention 
was paid to free-standing SLs. In addition, a number 
of calculations were performed for biaxially compressed 
(BaTi03)i(BaZr03)i and (BaTi03)2(BaZr03)2 SLs and 
for free-standing (BaTi03)„(BaZr03)„ SLs with n = 2, 
3, and 4. 

The calculations were performed using the density 
functional theory with pseudopotentials and wave func- 
tion expansion in plane waves as implemented in the 
ABINIT code.^*" As in the earlier study of these SLs^ 
the exchange-correlation interaction was described in 
the local density approximation (LDA). As pseudopo- 
tentials, the optimized separable nonlocal pseudopoten- 
tials constructed using the OPIUM program^^ were used; 
the local potential correction was added to them to im- 
prove the transferability. The parameters used for con- 
structing the pseudopotentials and other calculation de- 
tails are given in Refs. [3 and [iSJ. The cut-off energy 
for the plane waves was 30 Ha (816 eV); the integra- 
tion over the Brillouin zone was performed using the 
8x8x4 Monkhorst-Pack mesh for SLs with n = 1 and 
2 and the 8x8x2 one for SLs with n = 3-7. The re- 
laxation of the lattice parameters and atomic positions 
was performed until the Hellmann-Feynman forces be- 
came less than 5 • 10"^ Ha/Bohr (0.25meV/A). In the 
case of biaxially compressed SLs, the in-plane lattice pa- 
rameter was fixed equal to 0.97ao, 0.98ao, and 0.99ao, 
where uq = 7.7182 Bohr (4.0843 A) is the in-plane lattice 
parameter for the PA/mmm phase of the free-standing 
(BaTi03)i(BaZr03)i SL. The phonon spectra were cal- 
culated using the equations obtained from the density 
functional perturbation theory. The total spontaneous 
polarization in the SLs was calculated by the Berry phase 
method; its distribution over layers was calculated using 
the formula P^ — ^i^iZ* ^aUi^p from displacements Ui 
of atoms obtained in the polar phase and tensors of their 
effective Born charges ^*„fl in the paraelectric phase; 
here Wi are the weight factors equal to unity for atoms 
in the Ti-0 or Zr-0 layer under consideration, 1/2 for 
atoms in neighboring Ba-0 layers, and zero for other 
atoms. 
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FIG. 1. The phonon spectrum of BaZrOa in the cubic phase. 



The properties of BaTi03 obtained using the described 
approach were published previously^* and are in good 
agreement with experiment. 



III. RESULTS 



Structure of BaZrOs 



The structure of BaZr03 calculated using the de- 
scribed approach is in good agreement with experi- 
ment and results of earlier calculations i^^— For example, 
the lattice parameter in cubic BaZr03 is 7.8724 Bohr 
(4.1659 A) and differs from the experimental value 
(4.191 A at 2 K, Ref.0) by 0.60% (the small underesti- 
mate of the lattice parameters is characteristic of the used 
LDA approximation). As in earlier studies j^Ti^ the calcu- 
lations of the phonon spectrum of cubic BaZrOs reveal 
an unstable phonon mode at the R point at the bound- 
ary of the Brillouin zone with a frequency of 81i cm~^ 
(Fig. [T|). The unstable phonon at the M point observed 
in Ref. l9. is absent in our calculations. 

However, the instability of the phonon spectrum of 
BaZrOa at the R point has not yet been confirmed ex- 
perimentally. One of the possible causes of disagreement 
between experiment and theory can be quantum fluctu- 
ations^ which can destroy the long-range order in rota- 
tions of the oxygen octahedra.^^ According to our calcu- 
lations, the ground-state structure of BaZr03 is I A/ mem. 
This result differs from predictions of Ref. |^, in which 
the PI structure was considered to be equilibrium. In 
the present study, it was shown that the PI structure is 
nonequilibrium and slowly relaxes to the IA/m,cm struc- 
ture. The calculated static dielectric constant of bar- 
ium zirconate in the ground state is Sxx = £vv = 58.8, 
£^^ = 53.3 (the experimental value is 47, Ref. 
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FIG. 2. Phonon spectra of the free-standing 

(BaTi03)i(BaZr03)3 superlattice in (a) PA/mmm, (b) 
Amm2, and (c) Pmc2i phases. 



B. The BaTiOg/BaZrOg superlattices 

The phonon spectra of different phases of the free- 
standing (BaTi03)i(BaZr03)3 superlattice are shown in 
Fig. [2j One can see that two types of instabihty are 
observed in the phonon spectrum of the paraelectric 
PA/mmm phase (Fig. ^a): the ferroelectric instability 
and the structural instability associated with the oxy- 
gen octahedra rotation. The phonon frequencies at high- 
symmetry points of the Brillouin zone in the paraelectric 
phase of superlattices with n = 1, 2, and 3 are given in 
Table HI 

The instability region appearing as a band of imag- 
inary phonon frequencies along the path T—Z-R-X-T 
(imaginary frequencies are presented as negative num- 
bers in the figure) is caused by the ferroelectric insta- 



bility of ...-Ti-0-... chains, which was established for 
the first time Ref. [20. An analysis of the eigenvectors of 
phonons related to this instability region shows that out- 
of-phase transverse Ti and O atomic displacements in the 
xy plane in chains propagating in the [100] and [010] di- 
rections dominate in these eigenvectors at all points of 
the Brillouin zone; at the Brillouin zone center, this dis- 
placement pattern corresponds to the doubly degenerate 
ferroelectric Eu mode. For out-of-phase atomic displace- 
ments in the chains propagating in the [001] direction and 
consisting of alternating titanium and zirconium atoms 
(...-Ti-O-Zr-0-...), the lowest-energy ferroelectric ^2m 
mode at the F point is always stable; its frequency in- 
creases from 37 cm^^ for ti = 1 to 85 cm^^ for n = 3. 
A similar instability region along the path T-Z-R-X-T 
was observed in KNbOs/KTaOs superlatticesj^ 

Among the possible polar phases resulting from the 
condensation of the unstable E^ mode at the F point, the 
Amm2 phase with polarization along the [110] direction 
has the lowest energy in (BaTi03)i(BaZr03)„ SLs with 
n = 1-7. However, the phonon spectra of this phase 
(Fig. ^b) show that phonons at the M and A points 
at the boundary of the Brillouin zone remain unstable, 
which indicates the instability of the Amm2 structure 
with respect to the oxygen octahedra rotationi^i This 
instability of the polar Amm2 phases is characteristic 
of SLs with all periods (see Table H]). Thus, the Amm2 
phase is not a true ground-state structure for the SLs 
under consideration. 

As follows from Table U the frequencies of two unsta- 
ble phonons at the M and A points for the Amm2 phase 
of the superlattices with different periods are very close. 
Therefore, to search for the true ground state, it is nec- 
essary to compare the energies of all structures resulting 
from the condensation of these two phonons. 

An analysis of the eigenvectors of unstable phonons at 
the M point shows that the most unstable phonons are 
characterized by out-of-phase rotations of octahedra in 
neighboring layers. In SLs with even period {L = 2, 4), 
the rotation angle is small in the BaTiOs layers and large 
in the BaZrOa layers (Figs.Oa and[21c); in SLs with odd 
period (L = 3), the rotation angle is zero in the BaTiOs 
layers and is large in the BaZr03 layers (Fig. [31 6). Modes 
with lower instability (Table [I|) include modes with other 
combinations of rotations in layers (in particular, in- 
phase rotations in BaZr03 layers and strong out-of-phase 
rotations in every second BaZr03 layer) and a doubly de- 
generate mode with a complex combination of distortion 
and tilting of the oxygen octahedra (for L = 4). 

An analysis of the eigenvectors of unstable phonons at 
the A point shows that in SLs with even period {L — 2, 
4), the most unstable phonons arc characterized by the 
out-of-phase rotations in neighboring BaZrOs layers and 
the absence of rotations in the BaTiOs layers (Figs. [31rf 
and[3l/). In SLs with odd period (L — 3), out-of-phase 
rotation angles are large in the BaZrOs layers and small 
in the BaTi03 layers (Fig. [31 e). In both cases, the physi- 
cal period of SLs is doubled in comparison with the period 



TABLE I. The lowest frequencies of optical phonons at high-symmetry points of the Brillouin zone in PA/mmm, Amm2 phases, 
and in the ground state (Pm,c2i) of free-standing (BaTi03)i(BaZr03)n superlattices with n = 1, 2, and 3. 
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FIG. 3. The octahedra rotation angles in the ground state 
of (BaTi03)i(BaZr03)n superlattices corresponding to the 
condensation of phonons at the (a-c) M and (d-f) A points. 
{a,d) n = 1 (L = 2), (6,e) n = 2 (L = 3), (cj) n = 3 (L = 4). 
The BaTiOs and BaZr03 layers are denoted by open and filled 
squares, respectively. Vertical lines with arrows indicate the 
physical periods of the resulting structures. 



L specified by the layer sequence. Modes with lower in- 
stability (Table H]) are qualitatively similar to the above- 
described modes at the M point. 

An abrupt change in the rotation angle in going from 
the BaTiOa layer to the BaZrOa layer and its weak de- 
pendence on the period of SL correlate with the small 
change of the unstable phonon frequency along the M-A 
Hne, which indicates a strong localization of these rota- 
tions in the layers. 

To search for the ground state of SLs, the oxygen octa- 



TABLE II. The lattice parameters in the paraelectric 
PA/mmm phase of free-standing (BaTiOs )i(BaZr03)n su- 
perlattices with 1 < n < 7, the energy of the ferroelectric 
ground state A_Bi (the Amm2 phase), the energy A_E2 of the 
M/mcm or PA/mbm phases with pure rotations, the energy 
of the true ground state AE3 {Ima2 or Pmc2\ phase), and the 
polarization in the ferroelectric ground {Psi) and true ground 
(Pss) states. 
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4 
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7 


ao (Bohr) 


7.7182 


7.7747 


7.8009 


7.8160 


7.8258 


7.8379 


Co (Bohr) 


15.3298 


23.1923 


31.0604 


38.9304 46.8013 


62.5443 


AEi (meV) 


-72.7 
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-115.1 
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AE2 (meV) 
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-34.6 


-46.7 
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AE3 (meV) 


-78.3 
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-177.8 


-205.3 


Psi (C/m^) 


0.2755 


0.2024 


0.1568 


0.1273 


0.1070 
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Ps3 (C/m^) 


0.2591 


0.1859 


0.1445 


0.1113 


0.0984 


— 



hedra rotations corresponding to the least stable phonons 
at the M and A points of the Brillouin zone were added 
to the polar Amm2 structure. A comparison of the to- 
tal energies of the structures resulting from simultaneous 
condensation of the "rotational" and ferroelectric modes 
shows that the ground-state structure depends on the 
superlattice period L: for odd L, it is described by the 
space group Ima2 and results from the condensation of 
the unstable "rotational" mode at the A point; for even 
L, it is described by the space group Pmc2i and results 
from the condensation of the unstable "rotational" mode 
at the M point. The energy difference of the structures 
resulting from the condensation of phonons at the M and 
A points is rather small (1.6-1.9 meV). The energies of 
the structures corresponding to the ferroelectric ground 
state and to the true ground state, as well as to nonpo- 
lar structures with pure octahedra rotations (space group 



1.0 



0.1 



^ 0.01 



0.001 





















1 


, — ■ — 


n = \ 








/ 


N - A- 


« = 3 


^ 




A 


/ 


\ ••▼■■ 


« = 5 






_/' 




\S^ 








■ / 












/ 




























■.\ 








,▼ 




T^ 




r 


▼ 








y 




1 


1 1 




1 1 


, 







FIG. 4. The polarization distribution between the layers in 
the ground-state structure of (BaTi03)i(BaZr03)n superlat- 
tices with n = 1, 3, and 5. The layer number TV = corre- 
sponds to the BaTiOa layer. 



lA/mcm for odd L and PA/mbm for even L) are given 
in Table lllt^ An increase in all three energies with in- 
creasing n is due to two effects: (i) an increase in the 
BaZrOs volume fraction in the structure, which results 
in an increase of the instability of SLs with respect to 
the octahedra rotation; (ii) an increase in biaxial tension 
of the BaTiOa layers (as a consequence of an increase in 
the in-plane lattice parameter), which enhances the fer- 
roelectric instability. The corresponding change in the 
frequencies of unstable phonons at the F point can be 
seen in Table U 

The polarization profiles in the obtained ground states 
for SLs with n — 1,3, and 5 are shown in Fig. 21 One can 
see that the polarization is concentrated in the barium 
titanate layer and almost exponentially decreases with 
distance with a characteristic length scale of ^2.0 A. It 
is interesting that the spontaneous polarization in the 
true ground state is lower than the polarization in the 
ferroelectric ground state by only 6-12% (Table HIl). This 
indicates the weak effect of the octahedra rotations on 
the spontaneous polarization in SLs.— We explain this 
by the fact that regions with high polarization in the SLs 
under consideration are spatially separated from regions 
with strong octahedra rotations, and thus the compe- 
tition between the structural and ferroelectric instabili- 
ties becomes significantly decreased. This feature essen- 
tially distinguishes the superlattices from such crystals 
as CaTiOa and FbTiOs with competing instabilities in 
the cubic phase. In the former case, the structural dis- 
tortions completely suppress ferroelectricity; in the latter 
case, an opposite effect occurs. 

To test the possibility of forming the arrays of quasi- 
two-dimensional polarized planes in the ground state of 
(BaTi03)i(BaZr03)„ SLs, we estimated the interlayer 
interaction energy. The energy 2VFint (the energy cor- 
responding to two domain walls) was calculated as the 
difference of the total energies of ferroelectrically and an- 



tiferroelectrically ordered SLs with a doubled period in 
a similar way as in Ref. Ij; in this calculation, the octa- 
hedra rotations were neglected. The value of 2M^int was 
7.46 meV for n = 1, 1.505 meV for n = 2, and 0.420 meV 
for n = 3, and decreased almost exponentially with in- 
creasing n. For all n, this energy was much lower than 
the energy gain resulting from the ferroelectric ordering 
AEi (Table |TT| , which ensures almost independent po- 
larizations in the quasi-two-dimensional layers and the 
formation of the arrays of independent polarized planes 
similar to those observed in Ref. |j. The potential barrier 
height At/ for reorientation of the polarization, which 
occurs in the SLs under consideration by rotating the 
polarization vector in the layer plane, was 22.9 meV (per 
Ti atom) for n = 1, 33.2 meV for n = 2, and 38.3 meV for 
n = 3. The obtained potential barrier heights and ferro- 
electric ordering energies are much higher than those for 
the KNbOa/KTaOs SLs considered in Refs. 1 and i. In 
our opinion, these values are sufficiently large to provide 
stable polarization in the arrays of quasi-two-dimensional 
polarized planes at 300 K. 



IV. DISCUSSION 

We start the discussion with an analysis of possible 
polar structures in (BaTi03)i(BaZr03)„ SLs. The fer- 
roelectric instability in the ...-Ti-0-... chains can lead 
not only to the ferroelectric ordering, but also to the for- 
mation of antiferroelectric structures resulting from the 
condensation of unstable phonons at Z , R, and X points 
at the boundary of the Brillouin zone. Three types of the 
polarization ordering in the chains, which result from the 
condensation of these unstable phonons in the paraelec- 
tric PA/mmm phase, are shown in Fig. [SJ It is seen that 
these structures differ only in the relative orientation of 
polarization in neighboring chains. For example, for the 
phonon condensation at the X point, the antiferroelectric 
Pmma phase appears in the SLs; the modulation wave 
vector in this phase is directed along the x axis, whereas 
the out-of-phase displacements of Ti and O atoms form- 
ing the chains are directed along the y axis (Fig.[5l&)- It 
is remarkable that this phase still exhibits the ferroelec- 
tric instability with respect to the polar displacements 
along the x axis and finally transforms to the polar Pm 
{Pmc2i for n = 1) phase being ferrielectric. Similarly, 
the antiferroelectric Cmma phase with the modulation 
described by the phonon at the R point and with out- 
of-phase displacements along the y axis (Fig. [SJc) also 
is ferroelectrically unstable and finally transforms to the 
ferrielectric Abm2 phase. In contrast to non-degenerate 
unstable phonons at the X and R points, the unstable 
phonon at the Z point is doubly degenerate, and its con- 
densation can lead to two antiferroelectric phases, Pmma 
(Fig.[5l(i) and Cmcm, with displacements in one or both 
...-Ti-0-... chains, respectively. These displacements 
result in the in-plane microscopic polarization directed 
along [100] and [110]. It was the latter of these two phases 




FIG. 5. Patterns of the polarization ordering in the ...-Ti-O- 
... chains propagating along the y axis which accompany the 
condensation of unstable phonon at the (a) V, (h) X, (c) R, 
and (d) Z points. 



TABLE III. Energies E (in meV) of the antiferroelectric 
phases arising from the condensation of unstable phonons at 
the boundary of the Brillouin zone and of the ferroelectric 
Amm2 phase for free-standing (BaTi03)i(BaZr03)n super- 
lattices with n = 1, 2, and 3. The energy of the paraelectric 
PA/mmm is taken as the energy reference. 



Unstable 


n = 


1 


n = 


= 2 


n - 


-3 


phonon 


Phase 


E 


Phase 
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-37.3 


Pm,ma 


-52.1 
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-59.8 




Pmc2i 


-68.1 


Pm 
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Pm 


-109.0 
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Cmma 


-33.7 


Cmma 
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-59.7 




Abm2 


-64.9 


Abm2 
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Abm2 
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Pmma 


-41.9 


Pmma 


-60.7 


Pmma 


-69.6 




Cmcm 


-68.9 


Cmcm 


-100.3 


Cmcm 


-114.9 


V 


Amm2 


-72.7 


Amm2 


-101.0 


Amm2 


-115.1 



which was considered above in the calculation of the in- 
terlayer interaction energy. 

The ferroelectric instability of the intermediate Pmma 
and Cmma phases is easily understandable: in these two 
phases, the structural relaxation occurs in only one of two 
ferroelectrically unstable ...-Ti-0-... chains propagating 
in the [100] and [010] directions, and so the structure re- 
tains the ferroelectric instability in the second chain. In 
the Ptoc2i, Pm, Abm2, and Cmcm phases, the struc- 
tural relaxation occurs simultaneously in both chains; 
this is why these structures are stable (metastable). In all 
cases, the obtained intermediate and ferrielectric phases 
have higher energies in comparison with the ferroelectric 
Amm2 phase (Table [ml). 

Now we discuss the structure of rotations in the SLs 
under consideration. A comparison of the frequencies 
of unstable phonons corresponding to various combina- 
tions of rotations in the layers of SLs with the same pe- 
riod shows that phonons with out-of-phase rotations in 
neighboring layers are energetically least stable. These 
rotations are fully consistent with the instability at the 
R point in the parent cubic BaZrOa (Fig. [T]). The main 
difference between the two rotation systems correspond- 
ing to unstable phonons at the M and A points is the 
relation between rotations in the neighboring periods of 
the superlattice: for the phonon condensation at the AI 



point, when the structure is translated by one SL period, 
the octahedra are rotated in phase; in the case of the 
phonon condensation at the A point, they are rotated 
out of phase. It can be seen that the ground state of 
the SL always corresponds to such a rotation system in 
which the octahedra are rotated out of phase in any pair 
of neighboring layers, i.e., the M phonon is condensed in 
the SLs with even period and the A phonon is condensed 
in the SLs with odd period. 

One more subject for discussion is the comparison of 
the obtained results with the experimental data of Ref j-i^. 
According to the X-ray data from this work, the sponta- 
neous polarization in SLs is predominantly oriented in the 
layer plane; however, the line narrowing and an increase 
in the frequency of the soft E mode in Raman spectra in- 
dicate the monoclinic distortion of the structure and the 
appearance of the z-component of polarization in BaTiOs 
layers. The z-component of polarization determined from 
the hysteresis loops nonmonotonically depends on the su- 
perlattice period. The authors supposed that the appear- 
ance of the polarization component normal to the film is 
caused by the ferroelectric phase transition predicted in 
Ref. .6. in biaxially compressed BaZrOs layers. 

The calculations performed in the present work showed 
that the A2u phonon is indeed unstable in biaxially 
compressed BaZrOa (space group PA/mmm) with the 
in-plane lattice parameter equal to the lattice parame- 
ter of the free-standing (BaTi03)i(BaZr03)i SL (oq = 
7.7182 Bohr) (according to our data, the critical lattice 
parameter at which the instability of this phonon ap- 
pears is «7.735 Bohr and notably exceeds the value of 
7.425 Bohr reported in Ref. oj. However, the instability 
of the A2u phonon in biaxially stretched BaTiOa with 
the lattice parameter equal to ao is absent (the phonon 
frequency is 45 cm^^). Therefore, the appearing of the 
z-component of polarization in SLs needs a more detailed 
analysis. 

As we have shown previously,'^ the ^2u mode 
(which is responsible for the appearance of the 
z-component of polarization) is stable (37 cm~^) 
in the paraelectric PA/mmm phase of the free- 
standing (BaTi03)i(BaZr03)i SL. The calculations of 
the phonon frequencies for the PA/mmm phase of 
(BaTi03)„(BaZr03)„ SLs with n = 2, 3, and 4 per- 
formed in this work showed that the frequency of this 
mode rapidly decreases with increasing n (28.5 cm"-'^ for 
n = 2 and 17 cm~^ for n = 3) and becomes unstable 
at n = 4 (7i cm^^). However, the instability of the A2u 
phonon in the paraelectric phase does not mean that the 
z-componcnt of polarization will appear in the ferroelec- 
tric ground state. The in-plane spontaneous polariza- 
tion, which appears due to the strong instability of the 
Eu phonon, abruptly increases the frequency of the corre- 
sponding mode {Bi or B2) in the Amm2 phase, and this 
phase becomes stable with respect to the appearance of 
the z-component of polarization. 

This stability reduction of the A2U phonon with 
increasing n recalls the results obtained for the 



(BaTi03)m(BaO)„ Ruddlesden-Popper SLs^ In these 
SLs, a similar decrease in the A2U phonon frequency with 
increasing the BaTiOs layer thickness was observed; at 
TO > 8, it became unstable. However, already at m — A, 
the instability at the X point appearing in the SLs re- 
sulted in the formation of domain-like structure, since 
it was energetically more favorable to form a structure 
with a z-component of polarization periodically varying 
in space, rather than to uniformly polarize the BaO layer 
having a low dielectric constant. The macroscopic z- 
component of polarization is absent in such a structure. 

Another possible cause for the appearance of the z- 
component of polarization can be a substrate-induced 
strain in the superlattice (in Ref. il5i . the SLs were 
grown on MgO substrate). Calculations show that in 
(BaTi03)„(BaZr03)„ SLs with n = 1 and 2, for the in- 
plane lattice parameter equal to a = 0.99ao and 0.98ao 
(oo is the in-plane lattice parameter of the free-standing 
SL) and the disabled octahedra rotations, the polariza- 
tion vector is rotated by ~38° and ~67° from the layer 
plane (the space group of the ferroelectric ground state 
is Cm); for a = 0.97ao and disabled rotations, the vec- 
tor is rotated by 90° (the space group of the ferroelectric 
ground state is PAmm). However, when the octahedra 
rotations are enabled, in all SLs we see an unexpected re- 
sult: the z-component of polarization disappears in SLs 
with a — 0.99ao and 0.98ao (space group of the ground 
state is Pmc2i for n = 1 and Pnc2 for n — 2), whereas 
in the SL with a = 0.97ao, the polarization is rotated 
toward the layer plane by ^^24° (the space group of the 
ground state is Pc). 

The competition between the ferroelectric and struc- 
tural instabilities in crystals with the perovskite struc- 
ture has been known for a long time ;^^'^^ however, the 
fact that the octahedra rotations can have such a strong 
effect on the polarization seems to be observed for the 
first time. This finding can be important for the follow- 
ing reason. In the SLs under consideration, the energy 
gain resulting from the structural distortions is less than 
the energy gain resulting from the ferroelectric ordering 
(compare the energies AE'i and AE2 in Table IH)) . This 
leads us to expect that, as the temperature increases, 
the structural phase transition in the SLs will occur at 
a lower temperature than the ferroelectric phase tran- 
sition. Then, since the disappearance of the octahedra 
rotations results in the onset of the z-component of polar- 
ization, the situation is possible that an improper ferro- 
electric phase transition accompanied by the appearance 
of a nonzero z-component of polarization can occur in 
the SLs when the structural phase transition temperature 
is approached from below (when the octahedra rotation 
angle rapidly decreases). The abrupt rotation of the po- 
larization vector as the temperature varies can result in 
the appearance of a number of physical property anoma- 
lies resembling the anomaly in the piezoelectric coeffi- 
cients near the morphotropic phase boundary. However, 
currently, it is not clear whether the structural phase 
transition occurs in BaZrOa. Nevertheless, even if the 



establishment of the long-range order in the octahedra 
rotations in BaZrOa is impossible for some reason, the 
changes in the local angles of the rotations with temper- 
ature can have a significant effect on the z-component of 
polarization. 

The above-considered attempts to explain the ap- 
pearance of the z-component of polarization, which 
were based on the assumption of a uniformly strained 
BaTiOa/BaZrOa SL, were not successful (the strong bi- 
axial substrate-induced in-plane compression of the SL 
cannot be retained in a thick SL). It seems that the only 
way to explain this phenomenon is to consider the stress 
relaxations occurring in SLs with thick layers. This idea 
has already been used to explain the appearance of the z- 
component of polarization in BaTiOa/BaZrOa SLs with 
large periods. ■^•^'^ Our calculations show that the critical 
in-plane lattice parameter at which the z-component of 
polarization appears in BaTiOa is ^7.565 Bohr; in this 
case, the out-of-plane lattice parameter is c w 7.471 Bohr. 
If we take into account the systematic error (—0.7%) in 
the prediction of the lattice parameters for BaTiOa, the 
experimental parameter c — 3.981 A in the BaTiOa layer 
should correspond to this situation. According to the X- 
ray data,— such a relaxation of mechanical stress already 
appears in the SL with a period of L = 32 A. So, this 
relaxation explains the appearance of the z-component of 
polarization in SLs with the specified and larger periods. 

The first-principles calculations allow us to understand 
why the quality of short-period BaTiOa/BaZrOa and 
BaTiOa/Ba(Ti,Zr)Oa SLs grown in the [001] direction 
is not very goodJ^i^l This is because of the tendency 
of BaTii_a;ZrxTiOa solid solutions to three-dimensional 
chessboard-type ordering of cations at the B sites. The 
calculations show that the energy of the paraelectric 
phase of the (BaTiOa)i(BaZrOa)i SL with the elpaso- 
lite structure (space group Fm3m), which appears when 
growing the SL in the [111] direction, is by 83 meV (per 
formula unit) lower than the energy of the same SL but 
grown in the [001] direction. It is possible that defects 
generated during the stress relaxation in SLs grown in the 
[001] direction become the nucleation centers of three- 
dimensional chessboard-type ordered inclusions. Fur- 
thermore, the possible electrical activity of these defects 
can be a factor which changes the electrostatic boundary 
conditions at the heterointerface and stabilizes the polar 
state in the BaTiOa layer without inducing an apprecia- 
ble electric field and polarization in the BaZrOa layer. 



V. CONCLUSION 

The properties of BaTiOa/BaZrOa superlattices with 
competing ferroelectric and structural instabilities were 
calculated from first principles within the density func- 
tional theory. It was established that the quasi-two- 
dimensional ferroelectricity with the polarization ori- 
ented in the layer plane, which weakly interacts with the 
polarization in neighboring layers, occurs in SLs grown 
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in the [001] direction with one-unit-cell thick barium ti- 
tanate layer. It was shown that the energy gain from 
the ferroelectric ordering and the height of the potential 
barrier separating different orientational states of polar- 
ization are sufficiently large to observe the formation of 
an array of independent polarized planes at room tem- 
perature. The effect of the structural instability on the 
properties of SLs was considered. It was shown that 
the ground state is a result of simultaneous condensa- 
tion of the ferroelectric Fis mode and phonons at the 
M point (for SLs with even period) or at the A point 
(for SLs with odd period). Thus, the ground state is 
a structure with out-of-phase rotations in neighboring 
layers, in which highly polarized layers and layers with 
strong octahedra rotations are spatially separated. This 
significantly reduces the influence of structural distor- 



tions on the spontaneous polarization. It was shown that 
the switching on of the octahedra rotations in biaxially 
compressed BaTiOa/BaZrOa SLs with an in-plane strain 
of 1-3% results in an abrupt change in the polarization 
direction. This suggests that an improper ferroelectric 
phase transition accompanied by the appearance of a 
nonzero z-component of polarization can occur in these 
SLs with increasing temperature. 
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